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Abstract 

We investigate features of the lateral distribution function (LDF) of the radio signal emitted by cosmic ray air-showers with primary 
energies >0.1 EeV and its connection to air-shower parameters such as energy and shower maximum using CoREAS simulations 
made for the configuration of the Tunka-Rex antenna array. Taking into account all signihcant contributions to the total radio 
emission, such as by the geomagnetic effect, the charge excess, and the atmospheric refraction we parameterize the radio LDF 
This parameterization is two-dimensional and has several free parameters. The large number of free parameters is not suitable for 
experiments of sparse arrays operating at low SNR (signal-to-noise ratios). Thus, exploiting symmetries, we decrease the number of 
free parameters based on the shower geometry and reduce the LDF to a simple one-dimensional function. The remaining parameters 
can be fit with a small number of points, i.e. as few as the signal from three antennas above detection threshold. Finally, we present 
a method for the reconstruction of air-shower parameters, in particular, energy and Amax (shower maximum), which can be reached 
with a theoretical accuracy of better than 15% and 30 g/cm^, respectively. 
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1. Introduction 

The determination of the composition of the primary parti¬ 
cles is one of the most interesting and complicated problems of 
experimental high-energy cosmic ray physics. Imaging instru¬ 
ments, particularly, fluorescence or Cerenkov detectors, detect 
cosmic ray air showers with high precision, but their duty cy¬ 
cle is only in the order of 10%. On the other hand, detectors 
with a full duty cycle, such as particle detectors, until now have 
poor sensitivity to the shower maximum and cannot provide ac¬ 
curate studies of the composition. A candidate to solve this 
dilemma is the radio detection of cosmic rays. It probably can 
reach a precision comparable with air-Cerenkov measurements. 
However, it still has a number of important open issues such as 
efficiency, systematic uncertainties and precision of the energy 
and shower maximum reconstruction, all also depending on the 
detector layout. 

In the present paper we perform a detailed theoretical study 
based on a real large-scale detector layout. We performed about 
300 simulations based on the reconstruction of measured high 
energy Tunka-Rex Q and Tunka-133 El events. To simulate 
air showers we used CoREAS software integrated in COR- 
SIKA, which implements the end-point formalism for calculat¬ 
ing radio emission from air showers. In comparison with pre¬ 
vious work made for ideal detectors (see, for example, Ellll), 
LOPES 16] and LOEAR Qj our investigations have several im¬ 
portant differences. First, we reproduce detected events with 
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small uncertainty, thus, our simulation could be compared with 
signals measured by Tunka-Rex, which, in turn, features an ab¬ 
solute amplitude calibration. Second, the geometry of the de¬ 
tector matches modern large-scale setups, i.e. the spacing be¬ 
tween antennas is about 200 m. Finally, we transform the sim¬ 
ulated signals applying the real hardware properties of Tunka- 
Rex (amplifiers, antennas, etc.), and check the sensitivity of se¬ 
lected antennas. That means, we do a statistical study which 
gives realistic upper limits for the precision of the reconstruc¬ 
tion of air shower properties. “Upper limits” because we do not 
include noise and the precision will be slightly worse when tak¬ 
ing into account realistic background (see Appendix). There¬ 
fore, this limit could be reached in the case of large signal-to- 
noise ratios (SNR). 

The complication in describing the radio LDF originates 
from the interference of two completely different mechanisms 
of radio emission: emission due to geomagnetic deflection of 
charged particles, and the Askaryan (also known as charge- 
excess) effect. Adding these two effects causes an asymmetric 
two-dimensional lateral distribution function (LDF). There are 
two obvious approaches to describe this lateral distribution: to 
use a complex two-dimensional function, or to find some sym¬ 
metries and rewrite the LDF invariantly. The hrst approach was 
successfully tested in Q. It was shown that the LDF can be 
described with good accuracy, but the method used in this ap¬ 
proach requires a large number of points and, thus, a dense ar¬ 
ray. Because of that we used the second, customized approach. 
We found a transformation reducing the number of dimensions 
in the LDF representation to one, converting this function to an 
azimuthal-symmetric one. For this we estimate the strength of 
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the asymmetry based on the shower geometry. 

In additional to the simulation reproducing real Tunka-133 
events, we also performed an “ideal” simulation using a sym¬ 
metrical geometry and a three-dimensional dense detector. In 
contrast to the other simulations this is not for a statistical study 
of air-showers with different but for performing a tomo¬ 
graphic study of a mean air-shower (the description of this study 
is given in Chapter [2^ . In this way, we obtained evidence of a 
new feature of the radio emission, and found new connections 
of the parameters of the radio emission with the shower maxi¬ 
mum. 


Signal transformation and event selection on the detector 
level are made with the Auger Offline software framework il. 
We used the pattern of the Tunka-Rex antenna type, which is 
in first order close to a dipole. The frequency range is 30 - 
80 MHz. The event reconstruction pipeline is similar to Tunka- 
Rex, except for the SNR cuts: we do not add noise to the sim¬ 
ulation^ thus, we put only a threshold on the signal amplitude 
to reduce the digital noise. 

All plots, except Fig. are obtained with the Tunka-Rex 
layout and Tunka-133 event set. For Fig. [2|we performed a 


different simulation, as explained in Section 2.2 


1.1. Geomagnetic coordinate system 

To perform our calculations and later the reconstruction in 
an invariant way, we will use the so-called geomagnetic coordi¬ 
nate system, a special version of shower coordinates. The out¬ 
standing feature of this system is that the electrical field vector 
has only two non-zero projections to the axes, the third projec¬ 
tion is always close to zero. The basis of this coordinate system 


takes the form 

e.^VxB, (1) 

e, = V X (V X B), (2) 

e. = V, (3) 

where V and B are the shower axis and the Earth’s magnetic 
field (hat over vector means normalization, e.g. B = B/|B|), 


respectively. Let us also define useful angles: the geomagnetic 
angle - Z(V, B) and the geomagnetic azimuth 0g = /(e^, %). 

7.2. Simulation sets 

All simulations used in the present paper, have been pro¬ 
duced with CoREAS Q. As the hadronic interaction model we 
selected QGSJET-II. As the detector layout we used the setup 
of the Tunka-Rex experiment, which is located at an altitude 
of 675 m. The strength of the geomagnetic field was set to 
ai 60 pT, with inclination and declination of about 72° and 
-3°, respectively. Eor the incoming direction and energy we 
used measured Tunka-133 events from 2012/2013. We selected 
events satisfying the condition sin ag > 0.05 EeV. Tunka- 
133 reconstructs only air-showers with zenith angles < 50° 
due to design restrictions. That way, as initial parameters we 
used the energy of the primary particle £pr, the arrival direction 
(0, (p), and the core coordinates {x,y) on the detector plane. As 
the primary particle we used the two possible extreme cases for 
these energies: protons and iron nuclei. Due to the high resolu¬ 
tion of the Tunka-133 instrument, we can reproduce real events 
with high accuracy. The most important unknown parameter 
in the simulation is the depth of the shower maximum. Using 
different random seeds and primary particles (proton and iron) 
we try to limit the deviation between the shower maximum in 
the simulations and the measurement by Tunka-133 to less than 
30 g/cm^, as this is the precision of Tunka-133 12. Eor the 
present work we selected about 300 simulated events of each 
primary particle, using for each event the simulated shower with 
the smallest deviation between simulated and real X^ax- 


2. Asymmetry 

Presently there are a number of mechanisms for air-shower 
radio emission suggested by theorists Elll. We will consider 
only two contributions, which have been proven experimentally 
and which are the most important and dominant ones: geomag- 
netically induced transverse currents iTfOl and the Askaryan ef¬ 
fect ED. The complexity of adding these two contributions 
arises from the different mechanisms of the emission. While 
the electrical held of the geomagnetic emission is obtained by 
integration of charged particles Ne{h) over the height h and lies 
along the v x B vector, the Askaryan emission is mostly dehned 
by the derivative N'^ih) — dAe/dh and polarized along v - V, 
where v is the velocity of the particle and V is the shower axis. 
In Ref. fl^ good agreement between this simple model and 
measured data is shown. In our study we therefore assume that 
total polarization is a sum of two linear polarized contributions 
with unknown amplitudes. This leads to the known azimuthal 
asymmetry of the lateral distribution of the radio signal E2 El- 

While in principle the asymmetry can be extracted from po¬ 
larization measurements for individual events 02 El, in prac¬ 
tice, this is difflcult when the typical event has just three or for 
stations above the noise threshold. Therefore, we developed an 
approach to approximately correct for the asymmetry based on 
the event geometry only, since this is measurable with better 
accuracy for events with only few antenna stations. 

2.7. Origin of asymmetry 

The total electrical field at an antenna at each distance r and 
azimuth cpg can be represented as vector. 

£(r)=£g(r) + £ee(r)+£v(r), (4) 

where £g is a dominantly linearly polarized geomagnetic con¬ 
tribution, £ce is a radially polarized (like a normal Cerenkov) 
contribution from the Askaryan effect and £v ~ 0 is a vertical 
contribution to the signal. We neglect the contribution from the 
vertical component, since the angle between the shower plane 
and radio wavefront is only 1- 2° m- As in Refs. Bil, we 


^ This decision was motivated in order to avoid the study of the influence of 
noise on the signal, which is more appropriate in an analysis dealing with real 
measurements. Thus, the results we obtain here are a theoretical prediction for 
large SNR. For the influence of noise please see Appendix. 
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assume that the amplitude of £g and £ce changes only with dis¬ 
tance r, but is constant over (f)g. For £g also, the orientation is 
constant. Thus, the signal has the following components in the 
introduced geomagnetic coordinate system 

£g = (£g, 0, 0) = (£o sintTg, 0,0), (5) 

£ce = (£ceCOS0g, £ceSin(^g, 0), (6) 

where £g = £o sin ~ E^i- sin a^, and £ce ~ Fpi, and is the 
energy of the primary particle. We assume that the strength of 
the radio emission depends linearly on the energy of the elec¬ 
tromagnetic component (and, consequently, on the total energy) 
of the air-shower. The squared amplitude has the form 

£^ = (£o sin Ug + £ce cos <^g)^ -h sin^ (f>g 

- £q ((sin ag+ E cos (pgf' + sin^ (p ^, (7) 

where the asymmetry is dehned as e = £ce/£o|^ As we can see, 
when e > 0, the LDF is not azimuthally symmetric 




£>0 


Sir, (f>g) 


(9) 


To restore the azimuthal symmetry we dehne a transformation 
K (multiplication factor depending on r and pg) eliminating the 
azimuthal dependence. We can obtain the asymmetry e, if the 
two components of the electrical field are known 


Sx - So sin Og + £ce cos <f)g 
£y — £ce sin 0g 


( 10 ) 


This system can be solved for sin Ug sin 0g 0. Solving this 
system we obtain for the asymmetry e 

£ce £v/sin0g 

E - 


So Sx - £y cot 


■ sin ag 


( 11 ) 


Similar calculations have been also performed to extract the 
asymmetry from the experimental data of LOFAR ifTSll . From 
Eq.§ it is obvious, that K takes the form 


K = (, 
K£(r, 


£ + 


2s cos 0g sin ag + sin^ a. 


f 


g) ~ Scorr(0 “ So . 


( 12 ) 

(13) 


Although not used hereafter, as an aside, we want to note that 
the equation simplifies in the geomagnetic limit (£g » See), 
which was used as an approximation by several experimental 
analyses normalizing the measured amplitude by sin Ug 


K = 


1 


sinag 


» Scorr “ 


Sin Qg 


(14) 


Hereafter, we use the full form of K given in Eq. (121, not the 
approximation of Eq. (14 1 . 


^ In previous work Qinii similar calculations were performed with the 
notation 

a = sinagdCccl/lfigl) = £. (8) 

^ From the physical point of view, concerning the asymmetry, we discuss the 
interference of two components (geomagnetic and Askaryan). At ffg = 0 one 
of them disappears, consequently, term “asymmetry” has no sense anymore. At 
sini^g = 0 both components become collinear, i.e. £g = kSce, and cannot be 
distinguished at a single antenna. 


2.2. Asymmetry behavior 

Eor the asymmetry reconstruction we use the formulas given 
in Eq. ( [T0| ) to study the dependence of e(r) on the distance to 
shower axis r. The distribution of the asymmetry values for dif¬ 
ferent showers is broad. Thus we created a profile distribution 
for the average asymmetry values e(r) (see Eig.[^. The results 
for proton and iron are similar within about 10%. 



Figure 1: Averaged asymmetry profile of the radio lateral distribution for show¬ 
ers initiated by protons and iron nuclei based on about 300 CoREAS simula¬ 
tions of Tunka-133 events. Horizontal bars are the bin size, vertical bars are the 
standard deviation divided by the square root of the number of entries. 


Erom this picture we can conclude that the strength of the 
asymmetry varies not only from shower to shower, but also with 
distance to the shower axis. Similar results were obtained in 
Ref. ini. It was shown there that the asymmetry has a compli¬ 
cated structure, since geomagnetic and charge excess emission 
have their maxima at different altitudes, which leads to non¬ 
trivial addition of these components. Taking the longitudinal 
development of the air-shower into account, de Vries states that 
the behavior in the near region (up to 150 m from the shower 
axis) is mostly caused by Cerenkov-like effects. His statement 
does not explain the peak in the asymmetry appearing in this 
region. We give a potential explanation considering the lateral 
structure of an air-shower. The explanation could be the ex¬ 
istence of two sources for the charge excess contribution; in¬ 
tensive particle production and inelastic scattering close to the 
shower axis and particle absorption at distances far from the 
shower axis (the latter generally agrees with de Vries). In other 
words, in this peak we see the radio emission from the charge 
excess arising from very intensive particle interaction (like in 
the classical Askaryan effect in dense media). The behavior of 
the production and absorption regions could also depend on the 
distance to the shower maximum. 

To test this statement we studied a simulated radio prohle of 
an air-shower in detail: we simulated a vertical air-shower in¬ 
duced by a proton with energy lO'^ eV. The geomagnetic held 
was set to a strength of about 60 pT (similar to the strength at 
the Tunka valley) with geomagnetic angle o-g = 45°. The COR- 
SIKA simulated shower has its maximum at an atmospheric 
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depth of 636 g/cm^. We put several detector planes at differ¬ 
ent observation levels from 800 to 1000 g/cm^ with steps of 
10 g/cm^. Each layer consists of concentric rings with radii 
from 20 to 300 m with steps of 10 m. Each ring consists of 36 
antennas placed with azimuthal steps of 10°. By this we ob¬ 
tained a tomographic picture of the shower development in the 
region after X„ax- 

The obtained results show that the behavior of the asymme¬ 
try in the absorption region does not depend on the shower max¬ 
imum. But the position and height of the peak in the produc¬ 
tion region has a clear correlation with distance to the shower 
maximum (see Eig. |^. As we can see, it has the opposite be¬ 
havior than expected from the Cerenkov-like explanation given 
by de Vries El. In the case of a Cerenkov-like nature, the 
peak should move further away from the shower axis with in¬ 
creasing observation level, but we observe the opposite trend. 
Possibly, this indicates that after shower maximum the lateral 
extension of the shower region relevant for Askaryan emission 
shrinks closer to the shower axis. 

In the next step, we search for a simple parameterization 
of the asymmetry suitable for the reconstruction of measured 
air-showers. Although the production region depends on the 
shower geometry we neglect it, and approximate the asymme¬ 
try with a polynomial function ep(r) with fixed point £^(0) = 0. 
In the present work we do not pay attention to the relation be¬ 
tween the peak in the production region and the distance to the 
shower maximum. The peak in the production region depends 
on the shower geometry. Nevertheless, the relative height of the 
peak is small, especially for deep observers like Tunka-Rex (see 
Fig.|g, and, therefore, can be neglected in hrst order. The ab¬ 
sorption region has negligible dependence on shower geometry, 
therefore, we approximate the asymmetry as a function Sp{r) of 
distance to the shower axis r with fixed point £^(0) = 0 

£p(r) . (15) 

k>0 

Eor simplihcation we set = 0 and performed a global fit 
for both types of initial particles. The ht values are given in 
Table One can hnd the approximate mean asymmetry value 
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Figure 2: Top: Correlation between height and position of the asymmetry peak 
in the production region as function of the observation level of the shower. 
We marked the observation levels for several modem experiments for vertical 
showers: AERA, Tunka-Rex and LOFAR. The results are obtained for one ex¬ 
emplary vertical CoREAS shower with a primaiy energy of 0.1 EeV. The lines 
connect the points to guide the eye. Bottom: The distribution of asymmetry 
against distance to shower axis for two boundary observation levels: 800 g/cm“ 
and 1000 g/cm^. Highlighted points indicate peaks and show comespondence 
between top and bottom plots. 


Parameter 

Value 



(2.00 + 0.10)- 

10 ^ m ^ 


(-1.37 + 0.12) 

o 

1 

Ut 

3 

1 

hJ 


(4.10 + 0.46)- 

10-*^ m~3 

al 

(-3.67 + 0.56) 

- 10-" m-4 


Table 1: Global fit values for the asymmetry profile fit in Fig. □ of averaged 
CoREAS simulations. 

by solving the equation - Cn The asymmetry obtained in 
this way 3.1 r ^ 150 m is ^mean =^-11 ± 0.02. This value is in 


^ The statistical mean depends on the choice of the integration range, which 
can depend on shower geometry. To obtain a more stable value we decided to 
take the point of inflexion. This point is at distances of 100-200 m from the 
shower axis, where typically most antennas with signal are located, and where 
s is roughly constant. 


agreement with previous observations iHiia, where also the 
dependence of s over distance to shower axis and zenith was 
experimentally studied already ifTSll . It is important to note, 
that in spite of the geometrical invariance of s it stills depend 
on the strength of the magnetic field B, since £g = £g(B). That 
means, that the fraction of charge excess e should be roughly 
antiproportional with |B|. Consequently, the found value of e 
for the situation of Tunka-Rex is slightly smaller than the value 
found by the AERA and LOFAR experiments flAHTSll . 

The next question we studied was: what is the simplest 
function for description of the asymmetry s sufficient for a sat¬ 
isfactory description of the LDF after correction with K(£)? To 
test the quality of the correction we use a chi-square test. We 
define the goodness of the correction by the quantity e): 

the fraction of events passing the cut;if^/NDF < when the 
LDF is htted with K(£)£ 2 (^) (with &z{r) as defined in next sec¬ 
tion in Eq. ([T8]l). We start with the simplest function, a constant 
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( 18 ) 


value of the asymmetry ficonst- To find the optimal value for 
ficonst we solve the simple equation 

^jN^(Q^,e)dQ^^0 (16) 

The numerical solution of this equation gives the following val¬ 
ues: ^ 0.095 and ~ 0.075. If we compare the 

goodness of the correction made with this constant asymme¬ 
try and the parameterized one (see Fig. [^, we can see that the 
simple constant function gives a better result. From the one 
side, this can be explained by neglecting the peak close to the 
shower axis in the polynomial parameterization. The other im¬ 
portant factor is large deviations between individual simulated 
events at distances far from the shower axis (the spread at this 
distances is indicated by the larger error bars in Fig.[^, i.e. an 
effect of shower-to-shower fluctuations. A quantitative overes¬ 
timation of the asymmetry can even decrease the goodness of 
the correction. 

The result that a constant value of e yields on average an 
even better correction of the asymmetry than a polynomial makes 
the practical application of the correction approach very sim¬ 
ple. Consequently, we use a correction by a constant value 
s - 0.085 for the further analysis. 



N 

fr}(x) = ^ , 

k=l 


with parameters depending on the characteristics of the spe¬ 
cific air-shower. The parameter ro does not determine the shape 
of the function, and can be set to a defined, arbitrary value 
when fitting the LDF (we use ro = - 120 m for energy 

reconstruction and ro = re = 180 m for Xmax reconstruction). 
Using simple considerations (see Ref. ||3, Chapter 2.3.1) one 
can state that at the large distances the lateral distribution falls 
slowly, and the amplitude decreases below any detection thresh¬ 
old. This means, that the far-distance region is not useful for 
the reconstruction of the shower maximum. The slope func¬ 
tion would have no complicated features in the near-distance 
region a, if we did not take into account refraction in the at¬ 
mosphere (put refractive index rir = 1). But after introducing 
the refractive index we immediately obtain Cerenkov-like ef¬ 
fects imiiii. For the geometry defined in our simulations, the 
radius of the Cerenkov ring is at 100-150 m (see Fig. H}. This 
radius can be calculated from Eq. ( [T8] l solving the equation 

^£Ar(r) = 0, r = r(ro, Uk) (19) 

dr 

A solution can be found already with N - 2 



( 20 ) 


We make a comparison between different parameterizations and 
found that for the selected (Tunka-Rex) geometry a Gaussian 
(i.e N - 2) parameterization after asymmetry correction fits 
almost all event^ For N - 2, the goodness of fit N^(2^ = 
1) > 70%, while for a simple exponential LDF the goodness 
^x^Qx = 4) ^ 35% (see Fig.|gi. 

The properties of the Gaussian LDF such as mean (jj) and 
width (cr) are connected to the distance to the shower maximum 
(Fig.0 


cr = 


a\ 



1 


V-2a2 


( 21 ) 

( 22 ) 


Figure 3: Comparison between different methods of correction for about 300 
simulated showers. The fraction of accepted events N^(g^,f:) is calculated for 
the different forms of LDF K(e)£ 2 (^)- uncon'ected {s = 0), constant correction 
{s = Sconst = 0.085), correction with parameterization {s = £{r)). Using a 
constant value for the correction of the azimuthal asymmetry provides the best 
quality when fitting a one-dimensional LDF. 


3. Lateral distribution and its connection to the shower pa¬ 
rameters 

Already in the first observations of the radio emission from 


Now, that we found a good description for the lateral distri¬ 
bution valid for our detector, we tested methods for the recon¬ 
struction of air-shower parameters. We will follow the ideas de¬ 
veloped for optical air-Cerenkov emission, because after asym¬ 
metry correction the radio emission behaves similarly. 

The energy can be reconstructed by probing the signal am¬ 
plitude at a defined distance r^. Theoretical predictions for the 
optimal distance are about 50-150 m depending on the mass 
composition and geometry (see Ref. II 2 TII and later). Therefore, 
in general, the energy and LDF are connected by the following 
phenomenological relation 


air-showers an exponential-like function was suggested as LDF ifTSll . 
Later it was confirmed by modem digital experiments, e.g. LOPES lfT9l 
and CODALEMA If20l . Thus, we use an exponential-like func- 


F 

^pr 


= K 


/ fi(ro^re) r 
\ V/m j ’ 


(23) 


tion written in the form 

&N(r) = &ra sin ag exp[/;,(r - ro)], 


^ For fitting the LDF we used the true shower axis of the simulations as 
input. For Tunka-Rex measurement the shower axis is determined with only 
small uncertainty by Tunka-133. 
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Figure 4: Distribution of the quality of LDF fits for different parameterizations 
(Eq. |18) : exponential (A^ = 1) and Gaussian (A^ = 2). The fraction of events 
with LDFs corrected N^f^^/NDF, Sconst) and uncorrected N;^.{^^/NDF, 0) for 
the asymmetiy is averaged for both primaries (proton and iron). 
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Figure 5: Correlation between Gaussian LDF properties and distance to shower 
maximum. Average values for about 300 simulated showers. Definitions of ji 
(mu) and cr (sigma) are given in Eqs. \2i) and {22). 



where k is an amplitude slope parameter and is a power co¬ 
efficient. In other words, the energy is proportional to the am¬ 
plitude at a certain distance r^. To simplify experimental data 
analysis, one can set the arbitrary defined parameter rg in the fit 
to Te and take the fitted value £^0 with corresponding fitting un¬ 
certainty as energy estimator. To find we look simultaneously 
at the correlation between logarithms of energy and amplitude 
at different re and at the precision of the energy reconstruction 
using this formula (see Fig.|^. The maximum of the correla¬ 
tion points to the distance optimal for the energy reconstruc¬ 
tion. Since the relative difference between optimal distances 
re‘°'°‘’ and r“°‘\ and slopes and /c”"™ is about 10 % only, 

we selected median values for the energy reconstruction 

re = 120 m, (24) 


Figure 6: Comparison between true and reconstructed primary energy of the 
CoREAS simulations. 


EeV 


K = 422-, 

(25) 

V/m 

b = 0.93. 

(26) 


The precision of the energy reconstruction using these aver¬ 
aged parameters is better than 10 % for both particle types (see 
Fig.|^. 

For the reconstruction of the atmospheric depth of the shower 
maximum (Xmax) we use a slope parameter defined as 

d/„ &' 

77 = — = — = 2 a 2 (r - ro) -H oi. (27) 

dr o 
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Figure 7: Correlation between logarithms of amplitude and energy and preci¬ 
sion of the energy reconstruction using Eq.j^at distances from 40 to 160 m. 
The curve is the average of proton and iron simulations. 


rj is the slope parameter ai when evaluating the Gaussian LDF 
at the distance r = rp In other worlds, to obtain the slope rj of 
the lateral distribution, we take the htted value ai after setting 
ro - since the optimum distance for reconstruction 
is different from the distance re for energy reconstruction. In 
Refs, a nil a similar method was presented using the slope of 
the LDF. For reconstruction we use the parameterization 
suggested in Ref. Q 

^max = ^det/ COS 6 - (A + B log(ai H- b)). (28) 


This formula is more complicated than the one chosen for en¬ 
ergy reconstruction. It has two free parameters A and B which 
will be obtained from a fit to the simulated showers, one distance- 
dependent parameter ai = ri(r^), and one correction parameter 
b. We followed the same procedure as for the energy recon¬ 
struction: hnding the best correlation between the reconstructed 
and true shower maximum depending on the point {r^, b) in the 
two-dimensional space. One can see the correlation in a contour 
plot (Fig. 1^. 


As for the energy reconstruction, parameters in Eq. 27 have 
about 10% dependence on the particle type. After averaging the 
parameters and applying the formula, the relative difference be¬ 
tween true (simulated) and reconstructed Xmax values is smaller 
than 30 g/cm^ (see Fig.|^. The averaged parameters are 


Tx = 180 m, 

A = -1864 g/cm^ , B = -566 g/cm^ 
b = 0.005 m^' . 


(29) 

(30) 

(31) 



parameter & (xl0'^ m''') 


Figure 8: Contour plot of the correlation between LDF slope t] and true Xmax 
depending on the LDF parameter rx and the free parameter h. By fixing b, one 
obtains a distribution similar to Fig. Adding b as free parameter we obtain 
the two-dimensional dependence on {b,r^). In the chosen range of this two- 
parametric space, the correlation function behaves analytically and converges 
around {b = 0.003 m“\ = 195 m) and {b = 0.008 m“\ = 165 m) for 

iron and proton primaries, respectively. We have chosen the average point (b = 
0.005 m“^, Tx = 180 m) for further analysis, marked as star on the plot. 



Reconstructed distance to shower maximum (g/cm^) 


Figure 9: Comparison between true and reconstructed shower maximum of the 
CoREAS simulations. 


4. Conclusion 

We have found a way to describe the lateral distribution of 
the radio emission for a large-scale detector with antenna spac¬ 
ing in the order of 200 m. Our investigation of the azimuthal 
asymmetry has shown, that, in spite of the complex structure 
of interference between two different contributions to the radio 


emission (geomagnetic and Askaryan effects), a simple LDF is 
sufficient without introducing a large number of arbitrary pa¬ 
rameters. Assuming a constant value for the relative strength of 
the Askaryan effect provides a sufficient correction for the az¬ 
imuthal asymmetry of the lateral distribution. After correction 
for the asymmetry we need only a one-dimensional LDF with 
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3 parameters to describe Cerenkov-like effects. It was shown 
that for a non-dense radio detector array (with spacing in the 
order of hundred meters) using a Gaussian LDF is sufficient 
for primary energies in the EeV range. For the reconstruction 
of the energy and shower maximum we used formulas devel¬ 
oped for the air-Cerenkov detector Tunka-133. The free param¬ 
eters of these formulas were fit to CoREAS simulations made 
for Tunka-Rex events. The comparison between the true and 
reconstructed values of the simulations can be seen in Figs.|^ 
and|9] 

Our results are comparable with the precision of the Tunka- 
133 host experiment; the precision is better than 15% for the 
energy reconstruction and better than 30 g/cm^ for the shower 
maximum. The precision can be improved by applying qual¬ 
ity cuts. Let us note that about 20% of the uncertainty comes 
from the unknown chemical composition of the primary cos¬ 
mic rays (i.e. the precisions for a fixed particle types are about 
12% and 25 g/cm^ for energy and shower maximum, respec¬ 
tively). The reason could be a different distribution of Xmax 
for the proton and iron primaries, or different shower develop¬ 
ments. The developed methods will be used for the data anal¬ 
ysis of Tunka-Rex. The influence of noise has to be studied in 
more detail since the presented results do not include measure¬ 
ment uncertainties, e.g. due to noise. At high SNR the theoret¬ 
ically predicted precision should be achievable in practice (see 
Appendix). 

In future, our methods will be further optimized by exploit¬ 
ing information on the polarization, which is connected to the 
asymmetry. It is worth noticing that there is no evidence re¬ 
stricting our methods from being applied to more inclined events 
or events with higher energies, or other detectors with similar 
antenna spacing and frequency range. 
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Appendix: influence of background 


In the present Appendix we show the influence of the noise 
samples added to simulations on the reconstruction resolution. 


As input we used the dataset described in Section 1.2 and a 


library of noise measured by Tunka-Rex. The details of the 
pipeline used for reconstruction and quality cuts are described 
in II 22 I . in particular, we applied a cut on the SNR >10. 

The obtained results show, that noise has small influence on 
energy reconstruction after usual quality cuts (see Ref. Il22l ). 
However, these quality cuts have significant impact on the shower 


maximum reconstruction. First, almost all low-energy events 
are deselected decreasing the statistics of our dataset. Second, 
the resolution of the shower maximum reconstruction is slightly 
worse, about 40 g/cm^. This value is in reasonable agreement 
with measured data ll22ll . 



0.1 1 
Reconstructed energy (EeV) 


Figure 10: Comparison between true and reconstructed primary energy of the 
CoREAS simulations including realistic background. 
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Figure 11: Comparison between true and reconstructed shower maximum of 
the CoREAS simulations including realistic background. 
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